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I - Introduction 


This paper describes methods for determining spin-axis attitude 
(i.e., the direction in space of the spacecraft spin axis) and 
magnetometer biases which are being investigated for ground support of 
the Active Magnetospheric Particle Tracer Explorer (AMPTE) mission. 

The AMPTE mission will consist of two spacecraft.^ The first is the 
Ion Release Module (IRM), provided by the Federal Republic of Germany, 
which will be placed in a highly elliptical orbit with apogee at 
approximately 19 Earth radii in order to release lithium tracer ions 
outside the magnetosphere. This spacecraft will be spin stabilized at a 
rate of 30 rpm. The second spacecraft is the Charge Composition Explorer 
(CCE), which will detect the tracer ions inside the magnetosphere at 
altitudes of from 300 km to 7.5 Earth radii. The CCE will be spin 
stabilized at 10 rpm. 

Estimation of spin-axis attitude for both AMPTE spacecraft will be 
based on the measurements of the geomagnetic field and the projection of 
the Sun line on the spacecraft spin-axis, which we take nominally to be 
the symmetry axis of the spacecraft bus. 

For the purpose of this study, the attitude sensors are assumed to 
consist of a three-axis magnetometer and a Sun sensor which measures the' 
angle between the Sun line and For simplicity it is assumed 

likewise that one axis of the magnetometer is along The other 

two axes of the magnetometer define and I^. 

The measured quantities are taken to be 

M = magnetic field vector in body coordinates 

cos 0 = S*Y., where S is the unit vector directed from 

■~A 

the spacecraft to the Sun (6 is the "Sun angle"). 
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Attitude determination activities fall into two areas: 

• Determination of spin-axis attitude 

• Determination of the magnetometer biases 

Because the orbit -apogee distance for these two spacecraft is so 
great, accurate geomagnetic field data for attitude estimation is 
available only for the segment of the orbit near perigee. This is due to 
the poor accuracy of the magnetic-field model at such high altitudes 
resulting from both the small magnitude of the geomagnetic field as well 
as from fluctuations in the field caused by extraterrestrial phenomena. 
However, because of the large spacecraft angular momenta, it can be 
assumed for both spacecraft that the spin-axis attitude at apogee will 
not differ markedly from that at perigee of the same orbit. 

Algorithms for spin-axis attitude and magnetometer bias 
determination are now being investigated. These are: 

• attitude-independent estimation of three-axis 
magnetometer biases and 

estimation of spin-axis attitude from measurements 
of the Sun and geomagnetic field angle. 

Each of these algorithms are batch estimators utilizing a long segment of 
magnetometer and Sun data. The algorithms are developed in succeeding 
sections and then tested using simulated AMPTE data. 

II - Magnetometer Bias Determination 

The attitude of the spacecraft is usually not known before the 
magnetometer biases must be determined. Here an algorithm is developed 
which determines the magnetometer bias vector by minimizing a loss 
function which is independent of the attitude. 
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The quantities used throughout this section are defined as follows: 

H.(i) = jth component of the model magnetic field in the 
geocentric inertial (GCI) system at time i 

M.(i) = jth magnetometer reading at time i 

J 

For the ith point, an error <5(1) is defined by the following equation: 

^(i) = 


= jth component of the magnetometer bias vector, which 
is taken to be independent of the spacecraft 
position 


The objective of this equation is to minimize the quantity <5(i) by 
adjusting the bias vector^ to its optimal value. Thus, the loss 
function to be minimized is given by 


N 

L(B) = <^(i)|«(i) 


( 2 ) 


where w(i) is the weight associated with the ith data point. The weights 
are assumed to be normalized to unity, that is. 


I <^{i) =1 (3) 

i=l 

Determining the minimum value of L(B) first requires that its 
derivatives with respect to the components of the bias vector be set 
equal to zero: 


9L 

W~ = 


0 


m 


m=l ,2,3 


(4) 
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where 


Ci I 

L = ^ I u)(i) 


m 


i=l 


|H( 1 )|^ - |B - M(l)|' 


(B„-M„(i)) (5) 

m m' ' 


Combining Eqs. (3-5) leads to the following results: 


I G . B, = b + F (B) 
mk k m m 


(6a) 


or in matrix form. 


G B = b + ^(B) 


(6b) 


where 


B <Vk> 


(7a) 




= in <®m - V - 2 B-<M>B^ 


(7b) 

(7c) 


The bracket denotes the weighted average 


<A> = I w(i)A(i) 
i=l 


(B) 


6 , is the Kronecker delta defined as unity when m=k and zero 
mk 

otherwise. 


Eq. (6) can be solved directly to obtain the best value for the bias 
vector B, 


13-5 



General Description of the Iterative Solution 

Eq. (6) is nonlinear in JB and must be solved iteratively. The 
zero-th order (trial) solution to Eq. (6), is obtained by dropping the 
nonlinear terms in comparison to the linear terms. This approximation is 
valid only when the bias is small in comparison with the actual magnetic 
field. This point is not critical, as the iteration scheme constructs an 
accurate solution even when the trial solution is not close to the true 
solution. This will be discussed in more detail in the treatment of the 
numerical example. 

The trial solution is given by 

= G'^b (9) 

where G"^ = inverse of the matrix G 

= trial sol ution 

This solution may be iterated as 
The iteration continues until 



where e = some arbitrarily small value depending on the accuracy 
desired. 
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Numerical Examples 


The AMPTE engineering data simulator^ was used to generate biased 
magnetometer data for the purpose of investigating the convergence 
properties of the iterative solution. Two cases were considered: 

B/H « 1 


and 


B/H » 1 

The first case considered was B/H <1; in this case, 200 data points 

were used in the calculation. Data at the perigee point, at which the 

magnetic field attains its maximum value, was included. The magnetic 

field can be resolved into a component along the AMPTE spin axis, H|(, and 

a component perpendicular to the spin axis, Hj^. The maximum or perigee 

values for these components are H^^^ = 240 milligauss (mO) and 
MAX 

H|| = 90 mG. The input biases were chosen to be 5 mG, 10 mG, and 15 mG 

along the x, y, and z axes, respectively. The results of the bias 
determination calculation are shown in Table 1 taken from Reference 3. 
Rapid convergence and very high accuracy is obtained. The trial solution 
B^*^^ (iteration 0) initially was not accurate in the y component and 
needed to be iterated to obtain satisfactory results. Investigation of 
the case in which B » H used a subset of the data used in the first 
test. Here, 100 data points well outside the perigee region were used. 
For this test, = 5 mG and hJ*^^ = 2 mG. As before, the input biases 
are 5 mG, 10 mG, and 15 mG. These results^ are presented in Table 2. In 
this case, convergence is very slow and incomplete. Improved convergence 
cannot necessarily be obtained by using standard Newton-Raphson 
techniques. 
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ITERATION 

NUMBER 

LOSS 

FUNCTION 

(inG) 

By (mG) 

8^ (mG) 

0 

54621.0 

5.00288 

12.0278 

15.0213 

1 

5153.0 

4.98344 

9.38109 

14.9473 

2 

370.0 

5.00481 

10.1647 

15.0152 

3 

29.0 

4.99870 

9.95352 

14.9959 

4 

2.0 

5.00037 

10.0128 

15.0012 

5 

0.2 

4.99990 

9.99635 

14.9997 

6 

0.01 

5.00003 

10.0009 

15.0001 


Table 1 

Bias Determination Calculation for B/H < 1 


ITERATION 

NUMBER 

LOSS 

FUNCTION 

B (mG) 

X 

By (mG) 

B^ (mG) 

0 

24100.0 

1.8 

2.8 

5.3 

10 

1460.0 

3.7 

S.5 

11.0 

20 

501.0 

4.1 

6.1 

12.4 

30 

240.0 

4.4 

6.3 

13.1 

40 

133.0 

4.5 

6.5 

13.6 

50 

81.0 

4.6 

6.6 

13.9 


Table 2 

Bias Determination Calculation for B/H > 1 
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Ill - Spin-Axis Attitude Determination 


Once the magnetometer biases have been chosen properly, data from 
the Sun sensor and the magnetometers may be used to determine the 
spin-axis attitude. It is assumed that the spin axis is not varying over 
the data interval examined, 

A 

The spin axis is denoted by a^. The data are 


3(i) 

= measured Sun angle at time i 



= measured magnetic field at time i. 


l(i) 

= (true) Sun vector in GCI at time i, 
measured from the spacecraft to the sun 

I*”!)#** 

HO) 

= (true) geomagnetic field at time i , 

y N|y| 


Note that there will be no requirement of simultaneous Sun-sensor and 
magnetometer data. 

The spin-axis (attitude) vector, is subject to the following 
constraint: 

i-i»i (12) 

The spin-axis vector is chosen to minimize the following loss 
function: 


^S 


+ 


1 

7 


i=l 


<i). 


i( i ) 
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where 


X = Lagrange multiplier chosen to satisfy the constraint 
equation 

w^(i) = weight assigned to the jth magnetic field 
measurement 

oj|^(j) = weight assigned to the jth magnetic field 


measurement 


The quantity n is the angle between the geomagnetic field and the 
spacecraft spin axis given by 


n = cos-'(M /|h|) 


(14) 


The weights are normalized to unity 


Ns % 

i=l ^ i=l 


( 15 ) 


The spin-axis vector a is chosen to minimize the loss function 




( 16 ) 


m 


The derivative of the loss function is given by 


(i‘I(i) - cos e(i)) S^(i) 

m i = l 


(17) 


+ I <^M(i) (i*M(i) - cos n(i)) M^(i) - 


i=l 
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The solution to Eq. (16) may now be written as: 


,y''n,k - ' \ 


( 18 ) 


where 


= <S„S. >- + <M M. >„ 
mk m k S m k M 


= <cos S S„>- + <cos n 
m m S m M 


(19a) 

(19b) 


and the brackets denote weighted averages over the magnetometer and Sun 
data. That is. 


S 

<Cj>S "*s(i) Cj(i) (20) 

Eq. (18) may be written in matrix notation as 

(A - XI) a = b (21) 

where I is the unit matrix. 

Attitude Solution 

A general solution to Eqs. (18) and (19) is constructed in this 
section. The solution to these equations leads to the spin axis attitude 
in the Geocentric Inertial (GCI) coordinate system. Again an iterative 
procedure is developed to construct a numerical solution to the 
equations. An approximate solution to the problem is to take X = 0, 
i.e., to relax the constraint that a be normalized to unity. Given this 
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approximation, Eq. (18) may be solved to obtain 



( 22 ) 


Note that this vector is not normalized. In practice this solution will 
be very close to having unit norm since even with X = 0, a^ is overdeter- 
mihed in general by Eq. (18). Thus, normal i zing will lead to a 
very good approximation tor ^ (see Ref. 4). An exact numerical solution 
is generated by solving for X iteratively starting with a trial solution 
X = 0 and a^^^ given by Eq. (22). 

Define the function f(X) by 


f(X) = j(X)-a(X) - 1 (23) 

Given the numerical value of j(X), the Newton-Raphson method is used to 
determine X. Differentiating Eq. (23) gives 

l^=2a-^ (24a) 

and 


^ = (A - XI) ^ jL 

The Newton-Raphson scheme gives 


x(j) . 


f(x(j‘^^) 


(24b) 


(25a) 


a(j) = (A - x(j)l)-l^ 


(25b) 
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Numerical Example 


The spacecraft orbit in this example is of the AMPTE type, and the 
Sun and magnetometer data used covered the perigee point. The data is 
perfect (uncorrupted by random error) as generated by the AMPTE 
simulator. The "true" value of the right ascension, a, and declination, 
6 , were chosen to be 



ot = 159.67 deg 

(26a) 


<S = 0.0 deg 

(26b) 

The zero-order result as 

given by Eq. (22) was 



ot = 159.55 deg 

(27a) 


6 = 0.073 deg 

(27b) 

in very good agreement. 

After ten iterations. 

the values changed only 


slightly, as expected, namely 

ot = 159.76 deg (28a) 

6 = 0.062 deg (28b) 

IV - Conclusions 

Efficient and reliable algorithms have been developed for spin-axis 
attitude and magnetometer bias determination for the AMPTE spacecraft. 
Using simulated numerical data it was demonstrated that the methods work 
well for AMPTE mission parameters. The present work does not address 
problems associated with noise, data rate, sensor misalignments and etc. 
These problems were investigated in references (3) and (5). 
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